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Quantum simulations of High Energy Physics, and especially of gauge theories, is an emerging 
and exciting direction in quantum simulations. However, simulations of such theories, compared 
to simulations of condensed matter physics, must satisfy extra restrictions, such as local gauge 
and Lorentz invariance. In this paper we discuss these special requirements, and present a new 
method for quantum simulation of lattice gauge theories using ultracold atoms. This method allows 
to include local gauge invariance as a fundamental symmetry of the atomic Hamiltonian, arising 
from natural atomic interactions and conservation laws (and not as a property of a low energy 
sector). This allows us to implement elementary gauge invariant interactions for three lattice gauge 
theories: compact QED ([/(I)), SU{N) and Zat, which can be used to build quantum simulators 
in 1 + 1 dimensions. We also present a new loop method, which uses the elementary interactions 
as building blocks in the effective construction of quantum simulations for d + 1 dimensional lattice 
gauge theories (d > 1), without having to use Gauss's law as a constraint, as in previous proposals. 
We discuss in detail the quantum simulation of 2 + 1 dimensional compact QED and provide a 
numerical proof of principle. The simplicity of the already gauge invariant elementary interactions 
of this model suggests it may be useful for future experimental realizations. 
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I. INTRODUCTION 

The concept of quantum simulations goes back to the 
1980s, when Richard Feynman suggested the possibility 
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to simulate Quantum Mechanics using quantum com- 
puters pLj. Over the recent decades, the extraordinary 
progress both in theory and experiments, has enabled an 
unprecedent control in systems like cold atoms [11 13] , or 
trapped ions [4 . This control allows us nowadays to use 
them to simulate other quantum systems in an analog 
way, in the spirit of Feynman's visionary ideas. This has 
opened up a path to observe and understand many physi- 
cal phenomena which are, at least currently, unreachable 
in the context of analytic calculations or experimental 
measurements of the original systems [5 . 

Of specific interest are quantum simulations of many- 
body quantum models that appear in condensed matter 
physics. Some of those models are difficult to handle even 
with the most advanced numerical techniques, and thus 
quantum simulation appears as an important tool to in- 
vestigate them. These models include, for example, Hub- 
bard models, Heisenberg-like models in different lattice 
geometries. Atomic systems provide us with a natural 
playground to simulate those models, since once they are 
loaded in optical lattices, they are basically described in 
terms of simple Hubbard models whose properties can be 
tuned with external fields. During the last years, many 
quantum simulations of such models have been proposed, 
different techniques have been developed, and some of 
them have even been realized experimentally. 

Gauge theories constitute the fundamental building 
blocks of the standard model of High Energy Physics 
(HEP). They are built up out of fermionic and bosonic 
particles (or fields) , which represent matter and the force 
carriers, respectively. As a many-body quantum system, 
they are extremely rich in intriguing phenomena, and in 
some limits are very hard to study, even with the most 
advanced numerical techniques. Thus, a natural question 
is whether one could use the existing quantum simulators 
based on atoms in order to observe such phenomena, as 
well as to investigate regimes where standard techniques 
do not work. Unlike with condensed matter problems, 
however, the simulation of HEP models does not appear 
in a natural way in such atomic systems. In particular, 
they require: (a) Inclusion of both fermions and bosons; 
(b) interactions preserving local gauge invariance, which 
results in Gauss's law, as well as (c) Lorentz invariance 
(or its lattice remnants which ensure to proper contin- 
uum limit). 

In recent years, several works have proposed to use 
cold atoms in optical lattices to simulate HEP theories 
[6HT9]. Most of them propose to use optical lattices to 
enforce (c), and perturbation theory such that the desired 
interacting terms appear in the low energy sector. In 
some cases, this may lead to undesirable effects, since 
the gauge invariance is not exact, or since one has to go 
to high order perturbation theory, which leads to very 
weak effective interactions and strong constraints. 

In this paper, building on our early work on quan- 
tum simulation of HEP models [101 [131 [HI [HI [20] , we 
introduce new techniques in order to implement the con- 
ditions (a-c) above in an optimal manner. First, we pro- 



pose to use several internal states of the fermionic and 
bosonic atoms such that the gauge invariance of the re- 
sulting HEP model is a direct consequence of the con- 
servation of angular momentum in the original atomic 
scattering processes, and not a property of the low en- 
ergy sector after perturbation theory. Second, since the 
original Hamiltonian is already gauge invariant. Gauss's 
law does not have to be enforced: given that (in operator 
form) it commutes with the Hamiltonian, one just has to 
initialize the atoms in a state which satisfies it and then 
the dynamics will always occur in the subspace fulfilling 
that law. Thus, one could start out with the parame- 
ters corresponding to a regime where the ground state 
is well defined, and then turn adiabatically the physi- 
cal parameters in order to explore other regimes. Third, 
we propose to engineer the lattice system such that the 
traps for the bosonic atoms lie between the traps for the 
fermionic ones. In this way, fermionic hopping is medi- 
ated by a collision with the bosonic atoms, giving rise to 
the matter-gauge field interactions with the largest pos- 
sible coefficient, since the overlap integral responsible for 
this term is maximal. Fourth, we provide a new method 
(the "loop method") to construct the plaquette interac- 
tions that give rise to the dynamics of the gauge field (in 
the case of 2+1 dimensions). As in previous suggestions, 
we use perturbation theory in order to obtain the effec- 
tive terms in fourth order. However, we show that lower 
order terms in the perturbation series only renormalize 
our theory, so that we obtain the desired plaquette terms 
under the conditions that are equivalent to a second or- 
der perturbation (and not to fourth order). Furthermore, 
unlike in previous proposals, they are constructed out of 
already gauge invariant objects, and thus they do not 
require the explicit use of the Gauss law. 

The paper is organized as follows: in section II we put 
down the basics required for such quantum simulations: 
we brieffy discuss gauge theories and lattice gauge the- 
ories, and deduce the requirements that quantum simu- 
lations of such systems must fulfill; We discuss simpler 
high energy physics models which already show the inter- 
esting physics, but are simpler for quantum simulation, 
and review previous suggestions for quantum simulations 
of such systems. In section HI we describe the simulat- 
ing system - the general structure of the optical lattice 
and atomic Hamiltonian needed for such simulations, and 
then, in section IV, we show how to get, in the funda- 
mental Hamiltonian (without perturbation theory) the 
gauge invariant elementary interactions on links, for sev- 
eral gauge theories /7(1) , Z^r and SU{N). Then, in sec- 
tion V, we utilize these elementary interactions to build 
quantum simulations of 1 + 1 dimensional models - U{1) 
(a full nonperturbative simulation of the 1+1 dimensional 
Schwinger model) and SU{N). In section VI we intro- 
duce the new "loop method" and show how to use the 
gauge invariant nonperturbative elementary interactions 
to construct effectively plaquette interactions, required 
for the simulation of models in more than 1 + 1 dimen- 
sions, for these three gauge theories, and in section VII 
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we show how to use them and construct quantum sim- 
ulators for 2 + 1-d lattice gauge theories Zat and 
SU{N). The paper also contains an appendix, briefly ex- 
panding on some properties of the gauge theories whose 
quantum simulations are discussed. 



II. QUANTUM SIMULATION OF HIGH 
ENERGY PHYSICS 

A. Basics of High Energy Physics 

The standard model of High Energy Physics (HEP) is 
a Quantum Field Theory (QFT), in which the elementary 
particles can be divided into two separate groups. Matter 
particles (quarks and leptons) are fermions, represented 
by Dirac fields, while the force mediators, which are re- 
sponsible to the interactions among matter particles, are 
gauge bosons. Being a gauge boson, the gauge field must 
satisfy a special symmetry, which is Local Gauge Invari- 
ance. This symmetry may be either manifest or broken, 
but it is the nature of this symmetry which is responsible 
to the very special coupling between matter and gauge 
fields. Each gauge theory is based on a gauge group, 
whose elements are fundamental objects of the theory, 
forming the group space in which the gauge transforma- 
tions apply: these transformations do not correspond to 
changes in any physical observable, and thus they form a 
symmetry. The gauge groups may be either continuous 
or discrete, abelian and non-abelian. 

Both in abelian and non-abelian theories the gauge 
fields are massless [21]. However, in abelian theories the 
gauge fields are chargeless, while in non-abelian theories 
they carry charge. Quantum Electrodynamics (QED), 
for example, is an abelian {U{1)) gauge theory, whose 
charge is just the ordinary electric charge, and its gauge 
bosons, the photons, carry no charge. Quantum Chro- 
modynamis (QCD), describing the strong interactions, 
is a non-abelian (5/7(3)) gauge theory, whose charge 
is the color charge, carried by the quarks but also by 
the gauge bosons - gluons - due to the non-abelian na- 
ture of the theory. Abelian theories yield linear equa- 
tions of motion (Maxwell's equations for QED), since the 
chargeless gauge bosons do not interact among them- 
selves. In non-abelian gauge theories, there are such 
self-interactions, due to the non-abelian charge carried 
by the gauge bosons, and this results in nonlinear equa- 
tions of motion - the fundamental theory is described 
by the Yang-Mills equations for an SU{N) gauge the- 
ory [22 . Furthermore, non-abelian theories are respon- 
sible for long-range forces, manifested, for example, by 
the electromagnetic Coulomb law. However, non-abelian 
theories manifest the effect of confinement^ which binds 
matter particles together, such as quark confinement in 
QCD, which is responsible for the hadronic spectrum 
and forbids the existence of isolated free quarks [23l [24] . 
This non-perturbative phenomenon has been addressed 
over the recent decades in a variety of methods, includ- 
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FIG. 1. A part of the spatial lattice, in the k — 1 plane. The 
labeling of the vertices is shown. The links are labeled by 
their source vertex and their direction. For example, the link 
connecting between the vertices n and n + k is labeled as n, /c. 

ing Lattice Gauge Theories (LGT), where the space-time 
is discretized, enabling a numerical Monte-Carlo simula- 
tion [23, 25-27 . However, such a classical simulation, al- 
though very useful for many things - finding the hadronic 
spectrum, for example - still faces problems such as the 
sign-problem [28], which makes it hard to approach the 
limit of many fermions (a finite chemical potential), and 
thus probing some exotical phases of gauge theories (for 
example, color superconductivity in QCD) impossible 
using these methods [29H3T]. Other than that, quan- 
tum simulations enables also the simulation of dynamics, 
which is hard to simulate classically. 

Another important feature of HEP is being a relativs- 
tic theory, i.e. satisfying Lorentz invariance. This is 
of a great significance, as the interactions of elementary 
particles involve the regime of small distances and, of 
course, high energies, which requires a relativistic treat- 
ment, which is mostly avoidable in the case of condensed 
matter physics. While this symmetry must be exactly 
met in the continuum limit, it cannot hold in a discretized 
space-time as in LGTs. However, these theories still must 
include the proper remnants of Lorentz invariance, such 
that their continuum limit would be exactly relativistic. 



B. Basic ingredients of a Lattice Gauge Theory 

In lattice gauge theories, one can either discretize the 
entire (Euclidean) spacetime, or only the spatial direc- 
tions. As we are interested in a Hamiltonian model, we 
shall use the latter latticization, introduced by Kogut and 
Susskind [25H27]. 

In Hamiltonian LGTs, the fermionic (spinor) matter 
degrees of freedom, ip^^ reside on the vertices n G of 
a d-dimensional spatial lattice [32] (see figure [T]). These 
spinors may carry, generally, other indices, correspond- 
ing to possible physical quantum numbers of the matter 
fields, such as spin or flavor (which we avoid here) and 
also group space indices, which we denote here by lower- 
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case roman letters, a, 6, etc. And thus, generally 

4>n = (-i/'n.o) = I ^n,2 



(1) 



the dimension of the spinor depends on the representa- 
tion r of the gauge group used for it. Such fields may 
have local (mass) terms in the Hamiltonian, with the 
most general form 



M 



(2) 



where summation on the group indices is implicitly in- 
cluded, of course ('^j^'^n = ZlV^i a'^n,a )• Thcsc tcrms 

a 

are gauge invariant, as the group indices are fully con- 
tracted. Another way to see the gauge invariance is to 
consider the explicit local gauge transformation on the 
matter fields, 



b 



(3) 



where Vn is an element of the group, represented by the 
unitary matrix in the same representation of ip^^ de- 
fined locally for each vertex n, and see that it leaves these 
terms invariant. 

The interactions among the matter fields must include 
the gauge fields as well, being the force mediators. As 
such, the most reasonable place for the gauge degrees of 
freedom is on the lattice's links. Thus, on each link of 
the lattice, emanating from the vertex n in direction /c, 
define a group element Un,k (see figure [T|. These ele- 
ments can be represented by any representation r of the 
group. In general, Un,k are matrices of operators, de- 
fined on the link's local Hilbert space. This matrix space 
is called "group space", and the matrix indices, a, b etc., 
are referred to as "group indices". Explicit examples of 
such Un,k matrices of operators will be shortly presented. 

The interaction between neighboring vertices is me- 
diated using the link connecting them, in the form of 
elementary interactions (see figure [2| 



Hint (V^i^n,/cV^n+k + ^'C') 



n,/c 



h.c. 



n,k a,b 



(4) 



where is the matrix representation of U in r, the 
spinors' representation. Once again, these terms are 
gauge invariant as all the group's indices are contracted. 
One can also deduce, from the transformation law of the 
spinors ([3| and the gauge invariance demand, the trans- 
formation law of the [/'s, 



u: 



v:u: ,v 



n+k 



ab 



c,d 



(5) 



db 




n + k 



k 



FIG. 2. The elementary interactions - the interaction between 
the neighboring vertices n and n + k involves the gauge field 
on the link n, k connecting them: thus the gauge bosons are 
interaction mediators. As discussed in section [lV| this is the 
"natural" type of interactions in our simulation scheme, in- 
cluded in the fundamental atomic Hamiltonian. 



Finally, one shall introduce the pure-gauge terms as 
well. One type of pure-gauge terms which is gauge in- 
variant, is of the form 



Hp 



9_ 
2 



(6) 



where {E^^k)a generators of the group's alge- 

bra - for example, an angular momentum algebra for 
5/7(2), consisting of the angular momentum operators 
as generators. This term is just a sum of Casimir op- 
erators (which commute with all the generators, such as 
the total angular momentum operator for SU{2))^ and 
it is interpreted as the "electric energy". In general, 
one can define left {La) and right {Ra) generators on 
each link, constrained to give the same Casimir operator 

^LaLa = ^RaRa = ^EaEa. 

a a a 

These electric field also construct the generators of lo- 
cal gauge transformations, 

(Gn), = diVnK - Qn (7) 

where divn£^a is the "discrete divergence" of the group, 
div„E„ = ^((L„,fc)^-(i?„_fc,)J (8) 

k 

and Qn is the local charge (either dynamic or static). 
These generators are constant of motion - the physical 
states are the gauge invariant ones, satisfying, for each 
vertex n, the Gauss's law 



{Gn)a \phys) = 



(9) 



Other type of gauge invariant Hamiltonian terms is the 
trace of group elements along a closed path. The shortest 
such paths are plaquettes, and they form the "magnetic 
energy" part (see figure [3|, 

Hb = -\ Y1 (UiU2Ulul) + h.c}j (10) 

^ plaquettes 

where the 1,2,3,4 links are oriented along a plaquette 
(see figure [3|. The trace is on group (matrix) indices. 
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FIG. 3. The plaquette interactions - the gauge-gauge inter- 
actions. The labehng of the hnks around the pl aque tte is 
according to equation (10). As discussed in section 



VI 



these 

interactions are obtained effectively in our simulation scheme. 



Usually one includes in the Hamiltonian all such terms, 
where all the objects (group elements and spinors) are 
chosen to be in the fundamental representation. This 
will also be our choice throughout the paper. 

Next, let us give some explicit examples of three gauge 
theories we use and simulate in this paper. We shall 
describe the structure and Hilbert space of these theories, 
whereas further details can be found in the appendix. 



1. Compact QED 

As a basic example, we discuss compact QED (cQED). 
This is an abelian gauge theory, with the gauge group 
whose continuum limit is regular QED. However, 
unlike continuous QED, this theory manifests confine- 
ment of charges: at all values of the coupling constant g 
for 1 + 1 and 2 + 1 dimensions, and in the strong coupling 
regime for 3 + 1 dimensions [24l |26l |33l [34] . The com- 
pactness of the theory is essential for the existence of a 
confining phase [35 . 

In this theory, the /7n,/c operators defined on the links 
are pure phases: Un,k = e'^^^^^ . The conjugate electric 
field £^n,fe is merely an angular momentum operator, tak- 
ing integer values from — oo to oo. Thus, on each link 
Un,k the Hilbert space is the one of a quantum rotor, 
with canonical variables satisfying 



kl 



(11) 



this makes the U operators ladder operators of angular 
momentum - or, in other words, of electric flux: 



Un,k \m) 



m) 



1) 



(12) 



Note that as this group is abelian, there is no need to use 
different left and right generators. As there is only one 
generator. Gauss's law ([9| simplifies to 



(Gn) \phys) = 



(13) 



n+k 





FIG. 4. Illustration, using the t/(l), of the directionality of 
elementary interactions |15| If a fermion hops to the left, the 
flux in the middle increases. If a fermion hops to the right, 
the flux in the middle decreases. 



where Gvi 



n— k,/c 



Using these operators, we can deduce from the general 
Hamiltonians (|6|,([l0| the abelian version of the Kogut- 
Susskind Hamiltonian, 



Hrs — He 



Hp 



n,/c ^ n 



^n+i,2 



^n+2,1 



^n,2) 



(14) 



In the continuum limit. He is identified with the elec- 
tric energy and Hb with the magnetic one (as the 
cosine's argument is the curl of the vector potential: 

cos + 0n+i,2 " ^n+2,1 " ^n,2^ ^ ^ " ^) ' 

Hinti using staggered fermions [36 (see the appendix), 
we only have one spinor at each vertex, and in this case 
Q is simplified to 



Hint — e 



k 



(^/;te^*".'=V„+k + V'l+fce-'^-^^n) (15) 



and thus, the basic interaction involves a fermion hopping 
between neighboring vertices, while raising/ lowering, de- 
pending on the direction of the hopping fermion^ the elec- 
tric flux on the link connecting them (see figure [4|. 
For further details, refer to the appendix. 



2. Zn gauge theory 

Here we shall review the properties of a Hamiltonian 
Zn gauge theory [37]. We restrict ourselves to the pure 
gauge case, as only this is relevant for the purposes of 
this paper. 

First, let us describe the local Hilbert space on every 
link of the lattice. Define two operators, P, Q, which are 
unitary: 



Ptp = Q^Q = 1 



(16) 
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and satisfy the Z^v algebra, 



P^ = Q^ = 1 ; P^QP = e'^Q (17) 

where ^ = 

For example, one can work with the basis of P eigen- 
states. 



(18) 



with m e {-N/2, N/2} (WLOG, we assume N is odd 
- the change for an even N is straightforward) and then 
Q is a unitary ladder operator. 



Q\m) = |m-l) 



(19) 



with the cyclic property Q\-N/2) = \N/2). Altena- 
tively, one can expand the Hilbert space in terms of Q 
eigenstates, and then P will be a unitary raising operator 
(with the cyclic property, again). 

Interestingly, one can introduce the Hermitean opera- 
tors A on every link, by 



P = e'^^-Q = e' 



(20) 



and then, for ^ oo, one obtains the cQED Hilbert 
space, with canonically conjucate E,A. 

Let us now combine the entire lattice in order to get 
the gauge invariant Hamiltonian. It has the form 

H = He + Hb = 

n,/e 

- (Qn,lQ„+l,2QU,i<2 + h.c) (21) 

n 

One can define a static modular charge on the vertex 
n, by gn = e~^^'^ . Then the Gauss's law means, that a 
gauge invariant state must satisfy (for every n). 



Gn \phys) Qn \phys) 



where 



z+ z- 



(22) 



(23) 



with are links starting at n (positive links), and I— 
are ending there (negative links). 

Since these charges are modular and thus very different 
than the charges of continuous gauge theories, we shall 
only consider the pure-gauge case for Z^v in this paper. 



3. SU{N) gauge theories 

Let us first discuss the links' Hilbert space for a pure 
gauge theory. In a representation r, with representation 



matrices {T^}, the group elements can be parametrized 
as 



jjr 



(24) 



- it is a matrix in group space. 

Due to the non-abelian nature of the gauge group, it 
must have separate left and right generators, {La} , {Ra} 
respectively, corresponding to "left and right" non 
abelian electric fields [25, 21^- they can be represented as 
differential operators, canonically conjugate to the group 

parameters |^n,/c}- ^^^^ right generators of the 
group, they obey the following commutation relations 
with the group elements (within the same link, of course) 

[L„,C/n=r:[/'- ; [i?,,C/'-] = C/'^r: (25) 



and the Lie algebra 

[LajLb] = -ifabcLc ; [RajRb] = ^fabcRc 



(26) 



where fabc are the group's structure constants [38], and 
also [La^Rb] = 0. Physically, the difference between the 
left and right generators of a link may be interpreted as 
the color charge of it. The left and right generators can 
be obtained from each other using the group element on 
the link in the adjoint representation. 

From the local gauge transformation ([5|, one can con- 
clude that the generators of local groups transformation 
are ^ - where in the pure-gauge case Qn are C-numbers. 

From now on, we shall focus mostly on SU (2), the 
simplest continuous non-abelian group. There [25], r = j 
(total angular momentum quantum number), fabc = ^abc 
The local Hilbert space is characterized by three integer 
quantum numbers, j, m, m', which are eigenvalues of the 
Casimir operators and the z components of left and right 
angular momentum: 

^EaEa \jmm') ^LaLa \jmm') 

a a 

J2^aRa \jmm') = j {j + 1) \jmm') (27) 



Lz \ jmm') = m \ jmm') ; Rz \ jmm') = m' \ jmm') 

(28) 

the link Hilbert space may be interpreted as the one of a 
rigid rotator. The generators in the two edges of a link 
may then be interpreted as generators of rotations in the 
body/space systems [25^. The link's structure, in terms 
of operators and Young tableaux, is presented in figure 

El 

What shall be the Hamiltonian of such a theory? If 
we require it to be gauge- invariant, it may contain only 
gauge- invariant terms. Such terms can be constructed 
out of the generators, and since they must be contracted 
we get the Casimir operators LaLa = RaRa- They con- 
struct the local part of the Hamiltonian, called the Elec- 
tric Part - He (|6|. Other possibilities are closed loops: 
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n+k 



However, as we show, these terms can be obtained ef- 
fectively from the link terms, using perturbation theory 
[39l HO]. Although they are obtained effectively, gauge 
invariance is still fundamental, as the building blocks - 
elementary interactions, already fulfill the gauge symme- 
try. 



FIG. 5. The link's operators. Left and right generators 
{La} , {Ra}, and the group element. The generators trans- 
form in the adjoint representation (j = 1), as symbolized by 
the Young tableaux, and are related to each other by the ro- 
tation matrix U'^ (n, k) (which is in the adjoint representation 
as well). 



the "most local" ones are the traces of group elements 
directed around a single plaquette, forming the Magnetic 
Part We choose them to be in the fundamental repre- 
sentation (Ta = where a a are Pauli matrices) to 



get Hb (10), and finally obtain the SU (2) version of the 
Kogut-Susskind Hamiltonian [25^, H = He + Hb- 



Basic requirements for a HEP quantum 
simulation 



As can be understood from above, quantum simulation 
of HEP may be of great interest, and also significance, as 
it may help in avoiding problems of classical simulation, 
such as the sign problem. However, one should note it 
requires much more complex ingredients than quantum 
simulation of condensed matter systems - quantum sim- 
ulation of HEP models must: 

1. Include both fermions and bosons^ if one wishes to 
simulate both matter and gauge fields. This re- 
quires, for a cold-atom simulations, the use of many 
different atomic species. 

2. Respect local gauge invariance^ in order to have the 
correct symmetry which is responsible to the inter- 
actions and the interesting special features of the 
theories. 

3. Be relativistic. This can be reduced, if a lattice 
gauge theory is simulated, to the appropriate rem- 
nants of Lorentz invariance, but the continuum 
limit should still be relativistic. 

If one chooses to work on the lattice, as we do, the local 
gauge invariance "problem" transforms to the challenge 
of obtaining two types of interactions. First, the "link" 
gauge-matter interaction Q, which couples the matter 
and gauge field degrees of freedom in a very special way; 
Our basic idea is to get these interactions fundamentally 
in the atomic Hamiltonian - they will be derived directly 
from the conservation of hyperfine angular momentum F 
in atomic collisions. 

The second type of interactions is the plaquette in- 
teractions (10), which are, essentially, four-body interac- 



tions - not a fundamental part of the atomic Hamiltonian. 



D. HEP toy models 

Simulations of gauge theories, which must satisfy all 
the three requirements presented above, are challenging. 
However, when simulating HEP phenomena, one would 
not necessarily need, in first stage, to simulate the en- 
tire standard model, or even Quantum Chromodynam- 
ics. Several simpler models are available for observing 
the important phenomena and phases of the complicated 
theories. For example, working on the lattice. Compact 
QED is suitable for observing confinement (see appendix 
1): although everyday continuous QED manifests the op- 
posite behavior of a Coulomb phase, the compact lattice 
theory contains a confining phase in the strong coupling 
limit of the 3 + 1 dimensional theory, and confines for 
any value of the coupling constant in the 1 + 1 and 2 + 1 
dimensional theories [24, 26 , 33, 34 . Thus, for the obser- 
vation of confinement in a pure gauge theory, simulation 
of 2 + 1 cQED is enough (the 1 + 1 dimensional model is 
trivial). If one wishes to introduce dynamic charges, even 
the 1 + 1-d dimensional case is interesting - for example, 
one could simulate the lattice version of the Schwinger 
model [41 . 

As for non-abelian theories, full-fledged QCD with an 
SU (3) gauge symmetry is not essentially required as well, 
for the first step. A lot of theoretical, both qualitative 
and quantitative insight has been gained on QCD us- 
ing the 1 + 1-dimensional version of the theory, QCD2, 
or more generally, SU{N) in 1 + 1 dimensions [42H48] . 
On the other hand, some phenomena, such as confine- 
ment, may be observed also using a smaller gauge group 
- SU{2). Thus, for simulations of non-abelian gauge the- 
ories, SU (2) on the lattice [49|, even in 1 + 1 dimensions, 
is enough. 



E. Summary of previous works 

Several suggestions have been made for simulations of 
quantum field theories which do not include gauge fields. 
These include the observation of vacuum entanglement of 
a scalar field using trapped ions [50^, and the simulation 
of interacting scalar and fermionic fields - Thirring and 
Gross-Neveu models using cold atoms [6] (the latter could 
also be interpreted as a 1+1 simulation of fermions cou- 
pled to a gauge field). These two models correspond to 
simulation of fields in the continuum, respecting the ap- 
propriate relativistic and causal structure. Simulations 
of fermionic lattice QFTs have been proposed as well, 
where the fermions were either free or in external non- 
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dynamical gauge fields. These include Axions and Wil- 
son fermions [7 , Dirac fermions in curved spacetime jS] 
and general quantum simulators of QFTs and topological 
insulators [llj. 

As for abelian pure gauge theories, simulation of 2 + 
1-d cQED, with the possibility to observe confinement, 
first using Bose-Einstein condensates (BECs) of ultracold 
atoms in optical lattices [10 and then with single atoms 
in optical lattices [13 have been suggested, where the 
first is of the abelian Kogut-Susskind Hamiltonian [25] 
and the latter of a truncated "Spin-Gauge" theory. 

The inclusion of dynamical matter in such theories 
of great interest as well. This was done either for the 
link model |5TU54] - a 1 + 1-d simulation of the lattice 
Schwinger model [M], or as a generalization of the pre- 
vious pure-gauge simulations in 2 + 1-d, to include dy- 
namical fermions [TF. The latter also suggested a way 
to realize the gedanken experiment proposed in [20| of 
measuring Wilson-Loop's area law. 

All these abelian proposals fulfilled the relativity re- 
quirement through the use of the lattice. The models 
which included simulations fulfilled the first requirement 
by either including both fermions and bosons, or enabling 
the simulation of both types of particles. The gauge in- 
var iance demand has also been met, however, it has not 
been done in a direct way: gauge symmetry is not fun- 
damental in these models, but rather appears as a low- 
energy symmtery, manifested in the dynamics of an effec- 
tive Hamiltonian, obtained using a Gauss's law constraint 
required to introduce gauge invariance: this is since the 
four-body plaquette interactions are not fundamental for 
optical lattices. 

Simulations of other abelian lattice gauge theories are 
[m [55] . One should also note the continuum QED sim- 
ulation proposed in [9] (which does not manifest confine- 
ment as it is not compact [35 ). 

Some proposals for the simulations of non-abelian 
models have already been proposed as well, either uti- 
lizing prepotentials [56] , using ultracold atoms in optical 
lattices [T7| or utilizing Rishons in the link model [T8^ . In 
both methods a constraint is used in order to obtain the 
desired interactions. A digital simulation of an SU{2) 
pure gauge theory has also been suggested [19^. In [18] , 
as in non-abelian link models, the original symmetry is 
larger and one has to break it in order to obtain the 
right symmetry group; In [17 , the SU{2) gauge sym- 
metry is fundamental and is manifested already by the 
basic atomic Hamiltonian, unlike the effective methods of 
the former abelian simulations: this is done by exploit- 
ing the fundamental angular momentum conservation of 
the atoms, in a way which will be further explained and 
utilized for simulating other gauge theories in the next 
sections of this paper. 

A realization of discrete gauge theories (such as Z^) 
has been discussed using Josephson junctions [57j. 



F. The present work 

We have introduced the requirements from the quan- 
tum simulations of a gauge theory. As explained, in 
the previous proposals for simulations of abelian theo- 
ries, gauge invariance was effective, rather than exact. 
Here we shall describe the way to utilize a fundamental 
symmetry in systems of ultracold atoms in order to get 
a gauge symmetry which is not an effective low-energy 
symmetry, but rather built-into the theory, and thus is 
more robust. 

In the simulating scheme we suggest in this paper, we 
use fermions as matter and bosons as gauge fields, in 
vertices and links, respectively, like in the previous pro- 
posals. However, 

1. Local Gauge invariance. We do not impose Gauge 
invariance using an energy penalty in the Hamilto- 
nian. Instead, we show that by a judicious choice of 
fermionic and bosonic species (i.e., internal states), 
the natural atomic scattering interactions give rise 
to the terms we need with the appropriate gauge 
symmetries. This is so because the gauge symme- 
try in the resulting HEP model is equivalent to the 
angular momentum conservation in the collisions in 
the atomic model. 

2. Elementary interactions on links. The interaction 
terms between bosons and fermions are chosen so 
that they are maximal, and can compete with the 
real tunneling. This is obtained by using the idea 
of fig. [7|(see next section). 

3. Plaquette interactions. In 2+1 dimensional systems 
(and more), in order to obtain the dynamical terms 
of the gauge bosons (plaquette terms), we must use 
fourth order perturbation theory, introducing the 
loop method. This would naively mean that we get 
very small terms. However, we make sure that the 
odd orders are cancelled (or just renormalize pre- 
vious terms), so that in reality the conditions are 
equivalent to a second order perturbation theory, 
that are not so small. The plaquette interactions 
are O(e^) (where e is defined in eq. Q), but e, 
and the other odd orders of e are absent in the per- 
turbative series, and thus the expansion parameter 
is and effectively it is a second order contribu- 

tiom 0((e2)'). 
Resulting from that, we 

1. Propose a 1+1 dimensional cQED simulation, 
which should be relatively simple to implement ex- 
perimentally. 

2. Extend it to 2+1 dimensional cQED by adding pla- 
quette terms, and also introduce a new 2 + 1 dimen- 
sional model, Zat. 
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1. Single particle terms: 



FIG. 6. Schematic structure of the optical lattice used for 
simulations: Bosonic minima on the links (B), and fermionic 
minima on the vertices (F). 



3. Suggest a method for simulation of SU {N) theories, 
including a possible extension of the SU (2) model 
considered in 117 to 2 + 1 dimensions. 



III. THE SIMULATING SYSTEM 

Let us consider the atomic ingredients. We would like 
to build a theory of both fermions and bosons, where the 
fermions reside on the vertices, and the bosons - on the 
links (see figure [6|. Thus, let us start with the most gen- 
eral such structure. The vertices n of a square optical 
lattice coincide with the minima of fermions, described 
by the second-quantization operators (x) , where a la- 
bels the atomic species. Each link of this lattice coincides 
with a bosonic minimum, in which the bosons (x) may 
reside. If one assumes that the single-particle energy lev- 
els of each minimum are remote enough, only the lowest 
Bloch bands may be considered, and thus the second- 
quantized field operators may be expanded in terms of 
local annihilation operators Cn,^, an,/e,Q! and local Wan- 
nier functions '0n,a (x) , (/)n,fe,a (x), for fermions on the 
vertex n and bosons on the link emanating from it to 
the k direction respectively: 



(x) =^Cn,a^n,a (x) 

n,Q! 

(x) = ^ an,k,a(l)n,k,a (x) 
n,k,a 



(29) 



a ^ 

(30) 

where i^o,/, ^o,^ are the single particle Hamilto- 
nians, containing the kinetic energy and the trap- 
ping potentials, for the fermions and bosons respec- 
tively. Once the expansion (29) is plugged into 
these terms, and the overlap of Wannier functions 
is taken into account in the integration, one obtains 
two types of terms: local terms, linear in the atomic 
numbers, and nearest-neighbor hopping terms. In 
order to eliminate the latter for bosons, one should 
design the bosonic lattice deep enough such that 
any interactions outside a bosonic minimum would 
be negligible; In order to avoid fermionic tunnel- 
ing, one could use different species at neighboring 
vertices, alternately. 

2. Scattering terms: 

iisc = E sIm^ I (^) (^) ^7 (X) *5 (X) + 

a, (3,-^,5 

E a^^-ys J d'^^i (x) (x) *7 (x) ^5 (x) + 

a, (3,-^,5 

E 9a^,5 [ d'^H (x) (x) $t (x) (x) (31) 



Q!,/3,7,5 



where the scattering coefficients ^^J^^' 
are constrained by conservation laws and are fixed 
for different atoms - but can be controlled and 
modified using Feshbach resonances (perhaps op- 
tical, [5 8. -.60^ , if more than one is required). The 
first two terms represent the fermion-fermion and 
boson-boson scattering. Their integration, using 
(29), yields local scattering terms, within the same 



mmima. 



3. Rabi (laser) terms: 

HR=Y,^ip j rf'xvl/Kx)^'^ (x) 



a, 13 



Ef^f/3 / d'x$l (X)$^ (X) 

a, (3 



(32) 



The most general atomic Hamiltonian contains the fol- 
lowing terms: 



Using such terms, one can create "manually" de- 
sired hopping processes which may be useful in sev- 
eral cases. 

In principle one could also consider molecular terms, 
which dissociate into atoms, giving rise to a term with a 
bosonic operator and two fermionic ones. This may be 
useful for simulations in the bulk, but will not be used in 
the present paper where we concentrate on lattices. 
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FIG. 7. A schematic plot of the overlap of the fermionic Wan- 
nier functions (F) of two neighboring vertices and the bosonic 
Wannier functions (B) on the link. This is since the bosonic 
overlap is of order 1 and the fermionic tunneling is thus max- 
imal. 



IV. ELEMENTARY INTERACTIONS ALONG 
LINKS 

In this section we show how the fermion-gauge boson 
interaction terms appear in a natural way in the atomic 
system if one makes a judicious choice of internal states. 
Since the elementary interactions must come from the 
scattering of fermions with bosons (eq. (31)), it must 
involve an overlap integral between the initial and final 
bosonic and fermionic states. Given the fact that the 
fermions must hop, those two states will be located at 
different positions. In order to make this term as large 
as possible, one must have the bosonic atoms placed in 
between the fermionic ones (see figure [t]). Furthermore, 
in order to satisfy the gauge symmetry we will choose 
that the fermions and bosons change the internal states 
in this process according to the angular momentum con- 
servation. 

The key idea is angular momentum conservation: in 
these atomic scattering processes, the total hyperfine an- 
gular momentum Ftot is conserved. In particular, the z 
components - m^?, are conserved. One can specifically 
select the tuf values of the atomic species utilized, in 
order to generate the required interactions over the link, 
and eliminate the others. This will result in only gauge 
invariant terms, and forms the correspondence between 
two fundamental symmetries: angular momentum con- 
servation in the atomic, simulating level, is equivalent to 
gauge invariance in the simulated 

Let us first discuss the case of an abelian Hint, as in 



(15). 



A. U{1) elementary interactions 

For simulating cQED, we need two fermionic species 
and two bosonic species [61], arranged in an optical lat- 
tice, as in the previous section. Let us first consider a 
one dimensional lattice, and thus the links may be la- 



beled only by one index - the vertex from which they 
emanate - however, the same method may be generalized 
for more spatial dimensions, as we shall later do. 

We start with the bosons. Denote that bosonic species 
a, 6, both having two different value of m^. As explained 
before, no interactions take place between bosons of dif- 
ferent links. Thus the total number of bosons on each 
link is a constant of motion - we denote it by A/'o, setting 
it equal all around the lattice, and taking it to be an even 
number. 

On each link, a Schwinger algebra [62l [63] is con- 
structed from the two bosonic species: 



= a^b; L_ = b^a 



and 



N,);£=-{Na 



(33) 



(34) 



where Lz is our (truncated) electric field. 

Next, we wish to consider the fermions. As we would 
like to simulate a staggered fermions model ([36l[ll]5 see 
appendix 2), we only need a single fermion at most on 
each vertex. However, in order to use angular momentum 
conservation to ensure gauge invariance, we must use two 
different fermions, labeled by c and arranged such that 
the c minima occur in even vertices and the d minima 
in odd vertices. This eliminates the fermionic nearest- 
neighbor tunneling of Hq. The lattice is designed such 
that we get from Hq the mass Hamiltonian 



Hm = M^(-l)"^/;t 



(35) 



where t/jn is either Cn or d^, depending on the parity of 
the vertex. The Dirac sea state is obtained if initially 
all the d vertices are filled and the c vertices are empty. 
Note that if M > 0, the fermionic minima do not have 
to be mi?-dependent: if the system is initially prepared 
in a gauge invariant state, the fermionic tunneling of Hq 
is energetically forbidden and thus effectively eliminated 
and can be disregarded. Moreover, it also assures that 
two fermions can never occupy a single vertex (even of 
a different species), and thus the fermion- fermion scat- 
tering terms of Hgc may be disregarded. The fermions' 
local charges are defined as Qn = i^ni^n — \ (1 — (~1)"^) 
- for further details, see appendix 2. 

The gauge invariant elementary interactions are ob- 
tained from the boson-fermion scattering (the third part 
of Hsc ([31])). This is done by utilizing the total in 
atomic collisions. The hyperfine levels of the participat- 
ing atoms should satisfy 



rriF {a) -\- rriF {c) = tuf {b) -\- tuf {d) 



(36) 



(see figure |8| and thus, the only niF conserving processes 
(collisions) are (see figure |9|: 

1. a,c b^d^ and vice versa. This yields terms like 
clalbndn^i + dlbl^^an^iCn^2 
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ruF 





ArrtF 



FIG. 8. A schematic plot of the required choice of values. 
As in equation ( 36 ) , the equal spacing between Am f the m f 



of bosons (a, b on the right) and fermions (c, d on the left) is 
required to allow only "gauge invariant" collisions. 



(1) 




(2) 




FIG. 9. Schematic examples of the two types of possible 
Boson-Fermion scattering processes (collisions) along a link: 
(1) species changing and (2) non species changing. Both of 
them conserve total uif and thus are gauge invariant. 



2. a,c ^ a, c, and the same with b^d. This results 
in terms like cj^c^ (^'^k'^k + ^I^/e^ where the link k 
starts or ends in the vertex n. As all this terms 
a|.a/c + b\bk = A/'o, one eventually gets a contribu- 
tion which is proportional to the constant number 
of fermions - an ignorable constant in the energy. 

The first term is the desired gauge-invariant interac- 
tion. To see that, just perform the canonical transforma- 
tion 



and redefine the scattering coefficients in Hsc^ to obtain 



Hint — 



=^Y1 (V'i^+.nV'n+l + h.C.) 



(38) 



This Hamiltonian is especially interesting (although 
not realizable) in the limit Nq ^ oo. In that case, ^ ^ oo, 
and thus always m <^ i. Thus L±, in this limit, are uni- 
tary operators: 



\im) 



m (m ± 1) 



\£,m±l) 



\£,m±l) 
(39) 



and thus we get that ^^^^^^ approaches in this limit a 

unitary operator (pure phase), as in the Kogut-Susskind 
model. Another way to see it, is to consider that in this 
case the bosons form BECs. For TVq ^ 1, one can ap- 



proximate an 



vie 



etc. {m <^ i and 



thus Na — Nij <^ Nq and it is reasonable to approximate 



A^a^^,A^6^^Then, 



+1 



(40) 



where (j)n =6^ — 0^. This is similar to the mapping of 

m- 

Finally, we shall comment that the fermions can be 
replaced by hard core bosons on the vertices. 

B. Zn elementary interactions 

Now we turn to the construction of the elementary in- 
teractions of another abelian LOT theory, but this time 
with a discrete gauge group - Zn. Besides being an inter- 
esting gauge theory on its own, we consider its quantum 
simulation due to the fact that here, in order to simulate 
the exact theory, with exactly unitary gauge operators in 
the elementary intercations, we need a finite number of 
degrees of freedom, which makes this theory more tempt- 
ing to realize, unlike the cQED case, in which we only 
approximated the unitary interactions by angular mo- 
mentum ladder operators. 

Although we use the same general techniques of angu- 
lar momentum conservation, one must note that in this 
case it is not enough. This is due to the fa ct that the Zn 
Q operators are cyclic (see section II B 2| details), form- 
ing an "Escher's staircase" [65 , and thus regular angular 
momentum conservation is not sufficient. Therefore we 
use on top of the angular momentum conservation hy- 
bridization of states, and make use of auxiliary bosonic 
levels. 

For simplicity, we describe here the construction of el- 
ementary interactions of the Z3 case in 2 + 1 dimensions, 
but it can be easily generalized for larger TVs and higher 
dimensions. We do not consider the 1 + 1 dimensional 
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rriF 
— ► 



FIG. 10. A schematic plot of the required choice of val- 
ues, for the four "regular" levels {ai }^^ ^ , and the two "aux- 
iliary" levels {ci}^^2 equation ( |42| . Note that this plot 
shows TTiF only - not all the processes are available, only the 
"gauge invariant" ones, as in (45): some processes are elimi- 



nated due to energy shifts, which are not drawn in this figure. 
Also schematically shown (by arrows) are the lasers connect- 
ing coupled of levels in Hr ([46| 



theory as our fermions are not Z3 charges, but rather aux- 
iliary particles which shall be traced out in the derivation 
of plaquette interactions (see section VI B ): this point will 
become clear in the following derivation of the elemen- 
tary interactions. 

For obtaining the elementary interactions of Z3, we 
need, on each link, six fermionic species: four "regular" 
hyperfine levels, which we label {a^}^^^, and two "auxil- 

The vertices, unlike 
lose annihilation oper- 
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iary" levels {ci}-^2 (^^^ figure 
before, are occupied by bosons, w. 
ators are V^n,Xn which, due to energy shifts, can occupy 
alternating vertices (like the fermions in the Schwinger 
model simulation, and thus the use of "fermionic" let- 
ters). The vertex bosons are subject to a hard-core con- 
straint, 



He = A^TV, {N, - 1) 



(41) 



where Ny is the total number of bosons on the vertex 
V. Initially all the vertices are filled by only one boson - 
even vertices with ip and odd with x- 

The hyperfine levels of the atoms should satisfy the 
relation 

ruf (ip) + rrif (a^ = ruf (x) + ruf (a^+i) (42) 



for i G {1,2,3} (see figure 10). The c levels should be 
picked far enough energetically, such that they will not 
be involved in any link-species-changing process. Thus, 
the boson-fermion scattering terms will be of the the two 
following forms: 



1. Collisions with no change of species: 

{l,v) 



(43) 



where (l^v) are neighboring links and vertices, Ny 
is the total number of fermions on the vertex v and 
Ni - the total number of bosons on the link We 
prepare the system initially with Ni = 1. This is 
not changed by any interaction, and thus this term 
turns out to be proportional to the total number of 
fermions in the system - an ignorable constant. 

2. Collisions with a change of species. For even ver- 
tices (emanating from an even n), we have 



k \ i=l 



h.c 



(44) 



For odd links, one has to replace ^ x the equa- 
tion above. However, we can perform a canonical 
transformation, "inverting" the names of hyperfine 
levels on odd links, and then have them described 
by the same sort of interaction - compare to the 
canonical transformation of (37)). After doing that 
we call, formally, all the vertex bosons ip and obtain 



Hjr, 



n,/c 



V^iZ14n,/c«^,n,/c^n+k + ^.C (45) 



Since there is only one boson on each link, other boson- 
boson scattering processes are irrelevant. 



We also introduce, using Raman lasers (see figure 10), 



for each link, the following bosonic tunneling Hamilto- 
nian, within each link: 



Hr = (Ai +(5i) (^a\ai +aja4^+(l]i +a;i) (^a\a4 + ajai^ 
+ (A2 + S2) (a2a2 + 4<^2^ + (^2 + ^2) (cL2C2 + C2a2^ 
f (A3 + 53) (alas + 4^3) + (^3 + ^3) (4^3 + clas^ 



(46) 

Note that the ai ^ a^ process involves a three-unit an- 
gular momentum change, and thus it should be mediated 
by three photons. 

We make Ai^Qi the largest energy scales in the total 
Hamiltonian, and thus it will be reasonable to digaonalize 
Hji first, and obtain a hybridization of the couples of 
states coupled with lasers, in the form of a Bogolyubov 
transformation: 



(4 + 4); 
^(4 + 4); 
^(4 + 4); 



d 



't 



4 

71 ("^ 



(47) 



By setting f2j = — Aj, we obtain the diagonalized form 
Hr = 2Y, Ai4„,fc^i,n,fe + {5, Lo terms) (48) 



i,n,/e 
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and since we choose A to be very large, we can disre- 
gard, effectively, the di modes. Plugging the Bogolyubov 
transformation (47) into i^int^ disregarding the d modes, 
we get 



where 



Hint = (^i^n.fe^n+k + ^'^•) 



Q = b\b2^b%^blbi 



(49) 



(50) 



is the unitary Q of Z3 (see section II B 2 for details). Thus 
Hint is the desired Z3 elementary interaction. Note, how- 
ever, that the vertex bosons do not represent Z3 charges, 
and thus this method can only be used to generated aux- 
iliary particles, and not dynamic charges. 

Finally, we have to represent the electric part. He- 
Plugging the new modes into Hji^ we get 



'.,n,k 



(51) 



Then, setting Si = a;^, one gets 



(52) 



We identify 



P + pt = 2^ cos (mS) b^Jjm (53) 

m 

(where 5 = 27r/A/'), and thus, for Z3, 



(54) 

Set 5i = 5^ = ^ and 82 = — /i and obtain, neglecting 
constants. 



H = He + Hr + Hint = He + He + H^ 



int 



(55) 



- this is the fundamental Hamiltonian, with unitary ele- 
mentary interactions, from which we can now construct 
effectively the Z3 Hamiltonian with plaquette terms. 

Note that in order to obtain the correct interactions, 
one must use several Feshbach resonances. Their number 
can be reduced, if one generalizes the Hamiltonian to in- 
clude some energy difference between a[ai and a\a4^^ and 
also for the other two coupled of hybridized states. Then 
one can introduce a few more parameters to play with, 
and reduce the number of required Feshbach resonances. 

Also, note that in order to generalize to Zat for A/" > 3, 
one must have 2N bosonic species on each link, {a^}^^^ 
and {q}^2- The hybridization method is the same, 
with coupling between ai and aAr, and and q for 
iG{2,...,iV}. 
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FIG. 11. The lattice structure required for elementary SU{N) 
interactions. Each link is decomposed to left (L) and right (R) 
parts. It is constructed from two links of the optical lattice, 
tailored together by a constrained satisfied by auxiliary (aux) 
fermions in the middle. 



C. SU{N) elementary interactions 

We shall also refer to the fundamental Hamiltonian 
for SU{N) elementary interactions. There, the system 
is more complicated, and many atomic species are re- 
quired. Due to the decomposit ion of a single link to two 
parts (left and right, see section II B 3 ), this richer Hilbert 
space requires the construction of a single link of what 
was two separate links for the abelian theories - i.e., a 
simulating link is effective, and it is constructed from 
two atomic links, tailored by some constrained auxiliary 
fermions between them (see figure [TT]). 

We shall briefly review the ideas of [TT, in which a 
non-abelian quantum simulator for a 1 + 1 dimensional 
SU{2) lattice gauge theory was suggested. This simu- 
lation method utilizes prepotentials [56 , in which the 
group degrees of freedom are constructed out of "prepo- 
tentials" - harmonic oscillators, or, in our case, bosonic 
species. This enables a bosonic representation of the full 
Kogut-Susskind model. Fermionic representations are 
available too, using the link model [53l |54] - however, 
they correspond to truncated gauge theories, with finite 
local Hilbert spaces, from which one obtains the full the- 
ories only in the continuum limit. 

As explained in IIB 3[ and in figure [TTJ each link is de- 
composed into two parts, the left and the right, and hence 
simulated by two links of the optical lattice. In each of 
the link's parts, four bosonic species reside: ai,a2,ci,C2 
on the left, and 61, 62, c^i, c^2 on the right. The a, b species 
are the gauge field degrees of freedom, forming, using a 
Schwinger representation, the left and right generators of 
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the group, respectively, 

^« = l^^i ; = IYI^I {^a)kl k (56) 



satisfying the required algebra of the group (eq. 26), 



with j = 2 - 2 



= ^ and the Casimir operators lr = 



Nl (Nl 
2 \ 2 



Nr {Nr 



1) (where Nl 



a[ai 



b\bi + 62^2, satisfying the constraint 



Ul = 



-a2 
ai 



; Ur 



^ ( h\ 4 

-62 hi 



and obtain the group element on the link, in the funda- 
mental representation. 



U = UlUr 



(58) 



where 
Wl-- 



a[ci 
a\c2 



-a2C2 
aic\ 



Wr 



b\di 626^2 
—^2^2 hid\ 



(60) 



The auxiliary fermions are constrained by the large- 
scale energy constraint 



i),R^ 
a\a2 and Nr 
Nl = Nr). 

From these, in the prepotential method, one may con- 
struct the left and right matrices (of operators), in the 
fundamental representation. 



(57) 



satisfying the required commutation relations ([25|. 

The c, d species are prepared in coherent states (Bose- 
Einstein condensate) where a G M, a ^ 1. 

Let us denote the "real" spinors by ip and the auxil- 
iary ones, connecting between two links which will form 
one link in the simulated theory, x- Then, by prop- 
erly choosing the hyperfine levels of all the atoms (see 
the supplemental material of [17 for an explicit exam- 
ple), and tuning the scattering coefficients, we get the 
angular-momentum conserving interaction Hamiltonian 
for elementary interactions, 

+ ixl), .. + h.cA (59) 



and we label the two links from which the effective link 
n (emanating from the "real" vertex n) will be generated 
by n, L and n, R. 

The use of condensates for the auxiliary bosonic species 
allow us to replace Ci^di by a, and since a 1 we can 
approximately do the same for , , and one effectively 
obtains the Hamiltonian 

+ {Xi)i {UR,n),j {i'n+l)j VNR,n + l + h.C.] (61) 



>'^xiXn 



(62) 



If initially the system does not contain any x fermions, 
and A is the largest energy scale, we can obtain, using sec- 
ond order perturbation theory, an effective Hamiltonian, 
tailoring the two sides of each link, of the form 



h: 



eff 



V2 



+ 1 



h.c. 



(63) 

These are the elementary interactions of SU{2). How- 
ever, note that the bosonic link operators are not unitary, 
i.e. we have y/NL^n + ^Un\/ Nr^u + 1 rather than Un- In 
spite of that, as will be explained in the next section, one 
can still get qualitatively the same physics, in the appro- 
priate parameter regime. 

Also note, that although the full link is obtained effec- 
tively, the gauge invariance is still fundamental and it is 
constructed out of two already gauge- invariant building 
blocks: the left and right parts. 

The prepotentials method can of course be generalized 
to SU{N) gauge theories with N > 2 [66H68], and serve 
as a base for obtaining the elementary interactions in a 
similar manner. 



V. 1 + 1 DIMENSIONAL MODELS 

Having the elementary interactions in hand, we can 
now construct complete quantum simulations of 1 + 1 
dimensional gauge theories with dynamic fermions. This 
can be done for cQED^ but not for Z^r, as we do not 
discuss discrete charges in this paper. A proposal based 
on our method for the simulation of 1 + 1 dimensional 
SU{2) theory has already been suggested in y/Tj, and we 
shall review it here as well. 



A. Quantum Simulation of the Schwinger model 

Let us start with the a quantum simulation of 
Schwinger's model: a 1 + 1 dimensional abelian gauge 
theory (QED) coupled to dynamical fermions (see ap- 
pendix 2). Being 1 + 1 dimensional, this system does not 
involve any plaquette interactions, and thus we already 
have all the interactions we need. 

Besides Hint + Hm (eqs. ([^,([381)), we also need the 
electric Hamiltonian. He = = V5^^?,n5 

n,k n 

the atomic terms. 



Hk = 



-y 



{Nl^+Nl^-2Na,nN,,^) (64) 



This is exactly obtained from the boson-boson scattering 
terms of (31). These processes, of course, conserve the 
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total mp in collisions. The minus sign in the interaction 
may be avoided as well, thanks to the constant A^o* one 
could, instead, use the Hamiltonian 



2" - T ) 



(65) 



which is just He^ plus a constant in the energy - 
{Na^n + Nb,n) = <^Xl^o " which is, of course, ignor- 

n n 

able. 

Linear terms in the total number of bosons on a link 
(from Hsc and Hq) yield ignorable constants as well. 
Thus we get the Hamiltonian 



H = H 



E 



H 



M 



int 



(66) 



describing the dynamics of a U{1) "Spin-Gauge" theory 
with dynamic fermions [131 US] in 1 + 1-dimensions. 

For a finite A^o, one gets qualitatively the features of 
the model. As A^o (or ^) increases, the model becomes 
more accurate. The phase approximation can be made 
for condensates, in which one must make sure that three- 
body interactions are negligible. This can be assumed if 
the condensate is made in the shape of a tube, whose 
axis is perpendicular to the link, increasing the number 
of particles but reducing their density. 

Thus we have shown how to simulate a 1 + 1-d cQED 
with dynamic staggered fermions (lattice Schwinger 
model) using ultracold atoms, with an exact gauge sym- 
metry and no use of perturbation theory and effective 
low-energy considerations unlike in previous suggestions. 



B. Quantum simulation of 1 + 1 dimensional SU(2) 
gauge theory 

Having also the SU (2) elementary interactions in hand, 
one can obtain a 1 + 1 dimensional simulation of an SU (2) 
gauge theory. 



On top of the elementary interactions (|63|), one shall 
3lude as well t 



include as well the electric and matter Hamiltonians 



1 



2^V 2 V 2 
with gR + gL= 9^, and 



NR,n (Nr,„ , , 



(67) 
(68) 



This enables a simulation of the dynamics of the vac- 
uum of the theory, up to fifth order perturbation theory 
in H^^/ (63). See and its supplemental material for 
further details. 



VI. INTERACTIONS ON PLAQUETTES: THE 
LOOP METHOD 



In the next step, we would like to generalize our dis- 
cussion to further dimensions. However, the 1 -\- 1 — d ^ 
2 + 1 — (i transition is nontrivial, since the plaquette terms 
must be introduced, and, as explained, they are not a fun- 
damental part of the atomic Hamiltonian. In previous 
proposals, the plaquette terms have been obtained effec- 
tively, by constraining the Gauss's law, and introducing 
gauge invariance effectively, as a symmetry of the low- 
energy sector. In this section, we show yet another way 
to get the plaquette terms. Although in what we shall de- 
scribe the plaquettes will be obtained effectively as well, 
it is believed to be much more robust than the previous 
methods, since although we get the plaquette terms ef- 
fectively, gauge invariance is fundamental as described in 
the previous sections: the building blocks, which are el- 
ementary interactions, are already gauge invariant. This 
is called "the loop method". 

The idea is as follows. First, extend the scheme for 1 + 1 
simulations to more dimensions (which already serves as 
a simulation for the extreme strong limit). Then, treat 
the fermions as auxiliary particles, by adding a constraint 
which forces them to occupy only certain vertices (exactly 
one vertex belonging to each plaquette of the lattice). 
Hint, operating on states satisfying this constraint, will 
take us out of the "right" sector. Thus it would be rea- 
sonable to construct an effective theory in the ground sec- 
tor of this constrained Hamiltonian, and then, in fourth 
order (operating with elements of Hint around each pla- 
quette) one obtains the required interaction. Remark- 
ably, these fourth order terms are unnecassarily weak: as 
it turns out - and will be clarified throughout the follow- 
ing derivations - the relevant leading order is either the 
second order (for abelian theories, with no third order) 
or the fourth one (for non abelian theories), and thus the 
perturbative parameter should be small only to order 2 
or 4. 

The nature of auxiliary particles varies from one gauge 
theory to another, and depends on the gauge group. We 
shall describe, seperately, the methods of constructing 
such simulations for three different gauge theories: t/(l), 
Ztv and SU{N). For the sake of simplicity, we consider 
the 2 + 1-d case. However, at least by geometric means, 
the constructions for higher dimensions are similar. 

As a final general remark, before getting into spe- 
cific theories, one should note that in this method the 
fermions are "traced out" and eventually a pure-gauge 
theory is obtained. Such theories, in dimensions higher 
than 1 + 1, are interesting on their own, without includ- 
ing dynamic fermions; However, one could also include 
more fermionic species, not subjected to the "plaquette 
constraint" described above, which will serve as dynam- 
ical matter. We shall not elaborate on this possibility in 
this paper, but this is definitely possible and interesting. 
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A. cQED plaquettes 

Our first example of effective derivation of plaquettes 
in the loop method will be for the case of cQED - a gen- 
eralization of the Schwinger model simulation described 
in the previous section. For that, we start with a sim- 
ilar system to the one described for 1 + 1-d, but with 
two spatial dimensions instead of one. Thus, the system 
is described by the Hamiltonian ( 66 ) , generalized to two 
dimensions: i.e., all the vertex indices become vectors 
n — ^ n G and the links are now identified by two in- 
dices: n and the direction k G {1,2}. The electric part 
of the Hamiltonian is 



Yi.k 



(69) 



We introduce another set of fermions, Xn, which be- 
have like the ?/^nS (including the same interactions with 
the bosons). At this stage, for the sake of illustration of 
the method for obtaining plaquettes, we assume that the 
bosonic link operators are really unitary {Nq oo), i.e. 
we work with the elementary interactions 



Hint = (V^i^n,/c^n+k + Xn^n,/cXn+k + ^'^') = 

n,k 

eE (^^.e'^-'^n+k + xJ.e'*-'=X„+k + h-c.) (70) 

Note that the hyperfine levels of the two x species must 
be chosen carefully, such that X ~ remain separate in 
Hint and no mixing interactions can occur. The differ- 
ence between the types of fermions is found in a con- 
straint we add on the fermions. 



■ Fx (n) xlXn) 



(71) 



where is zero everywhere, unless where both the in- 
dices ni, n2 are even, where it takes the value of 1, and 
is zero everywhere, unless where both the indices ni,n2 
are odd, where it is 1. If we define both these types of 
vertices as even ones, we see that He puts an "energy 
penalty" for each species not being in its specific "pre- 
ferred" type of an even vertex (see figure 12). 

We denote the ground sector of He as A^o, and wish 
to work in this subspace. Thus, the system has to be ini- 
tially prepared in a state where all the fermions occupy 
only even vertices (the opposite to the Dirac sea case for 
dynamic fermions). Note that since H is gauge invariant, 
and we initially prepare the system in a gauge invariant 
state, the dynamics will leave the state gauge invariant 
and thus we choose to specifically work in M C Mo^ 
which is the set of gauge invariant state inside Aif. As A 
is the largest energy scale, we derive an effective Hamil- 
tonian within A^o - we shall construct a low-energy effec- 
tive theory which includes the plaquette interactions. In 
order to do that we use time-independent perturbation 
theory, following the notations of [39j. 



X 



DC 



FIG. 12. The "preferred" vertices of the V^, x fermions, in the 
even vertices, according to the constraint set by He (equation 
(71)). Each plaquette contains exactly such two vertices, one 
of '0 and one of x- 



1. First and second order contributions 

Denote Vq as the projection operator to A^o, and de- 
fine Hi = He ^ Hint and 



JC 



E 



\a) {a\ 



Ec {a) - Ec (0) 



(72) 



where Ec is the eigenvalue of Hc^ and thus Ec (0) = 
0. Also, for the convenience of series expansions, denote 



H^ 



The first order term in the effective expansion is 
VqHiVo = He- In second order, we have 



(1) 

eff 
(2) 
eff 



-VqHiJCHiVq. Here, in /C, only Hi will con- 
tribute, taking to (and from) intermediate states \a) 
with Ec (o^) = A (the constraint is violated only for one 
fermion). The contributions will be only of double oper- 
ations of Hi on the same link (see figure [l3^) , and due 
to the unitarity of the interactions will lead(within Mo) 
to a constant in the energy, which is ignorable. 



2. Third order contributions 

The third order contribution takes the form H^^j = 

VoHiJCHilCHiVo - ^{VoHiJC^HiVo^VoHiVo}. The 
second (anti-commutator) term is just a combination of 
the first and second order terms, which will result in 

2 

—Mj2He^ where J\f is the number of links. 

The first term is nonzero only for the combination 
VoHintfCHEVoJCHintVo' Here, as in the second order, 
Ec (o^) = A. These terms will vanish, unless we consider, 
as in the second order, double operation of Hint on the 
same link. Define a "positive" link if it starts on an even 
vertex, and a "negative" link if it ends there. Since only 
states where even vertices are occupied belong to Afo, 
only the part V^^^j^^n fc^n of Hint acting on Afo will 
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Hamiltonian, of the form 



FIG. 13. Examples for "jumping" of auxiliary fermions from 
their "preferred" vertices and back, (a) in the second order of 
Hint, (b,c) in the fourth order of Hint, where (c) forms the 
plaquette interactions. 



give rise to a nonzero contribution, and thus in the final 
operation, the contribution will come from V^i/7n,/cV^n+k 
(and similarly for x)- For negative links, only the op- 
posite processes contribute. Thus, for each positive link 
n, /c we get the contribution 



A2 



UHeU^ = 



(73) 



where C is the same constant as in the other third order 
contribution, and V is another constant. Similarly, for 
negative links we get 



En,k + V 



(74) 



The V parts can be ignored, as constants, and after 
summing on all the links (remember that each link is a 
neighbor of exactly one even vertex, and thus it is either 

2 

positive or negative) we get a cancellation of the j2He 
parts, and eventually we are left, for every even n, with 
^^divn£^, where the discrete divergence is divn£^ = 
[e^^^i + £;n,2 - ^n-i,i ~ ^n-2,2) const, within Mq 
by Gauss's law (in case of no dynamic charges). Thus 
all the third order contributions are constants and we 
disregard them. 



3. Fourth order contributions 



A3 



(75) 



n,/c 



- this is the reason for the inclusion of x fermions; Oth- 
erwise, we would have got such a renormalizations only 
for links which are neighbors of the "preferred" vertices 
of ips. 

The other type contributions are the anticipated pla- 
quette interactions. These arise from the contribu- 
tion -VoHintJCHintJCHintJCHintVo, where the four oper- 
ations of Hint cause a fermion to hop from its "rest" ver- 
tex, around a plaquette, and back to the original place, 
and thus forming the desired Hamiltonian terms (see fig- 
ure [iS]^). Each plaquette operator is divided into two 
different orientations (clockwise and counterclockwise). 
Each orientation is obtained twice, because it has to start 
(and finish) in an even vertex, and each plaquette con- 
tains two such vertices - i.e., each plaquette interaction 
is obtained once using ips and once using xs. Altogether 
we get the required 

n 

COS + <?^n+i,2 " ^n+2,1 " ^n,2 j (76) 



4e^ 
A3 



and eventually the abelian Kogut-Susskind Hamiltonian 
is obtained 



Heff = He + SHe + He 



(77) 



One may claim that the plaquette interactions are neg- 
ligible here, since they are obtained in fourth order per- 
turbation theory. However, note that all the terms in 
the effective Hamiltonian involving e actually involve - 
only even orders contribute (since one can return to Mo 
only with an even number of Hint operations). Thus, it 
is sufficient to demand <C A^, rather than e <C A. 

Dynamical fermions may be introduced through the 
inclusion of another set of fermions, with its own Hint 
and i^M, and without any constraint. 

In section IVIII we discuss a similar simulation in the 
real case - i.e., under real conditions, without the ideal 
Nq ^ 00 assumption, along with a numerical proof of 
principle. 



In fourth order, there are much more contibutions. The 
ones involving He yield mostly ignorable constants based 
on Gauss's law, as in the third order. There are some 
possible fourth order processes involving only Hint on 
two links, corresponding to back-and-forth hopping of 
fermions (see figure 13 3), which result in constants due 



to the unitarity of interactions, as in the second order. 

The non-constant contributions are of two types. The 
terms involving He yield a renormalization to the electric 



B. Zat plaquettes 



The effective construction of the Tjm Hamiltonian out 



of the fundamental Hamiltonian (55) in the loop method 



is similar to the derivation in the cQED case. The anal- 
ogy applies to the ideal limit of the cQED simulation, 
as the interactions here are exactly unitary (i.e., the link 
operators are unitary). As before, one sets A to be the 
largest energy scale, comparing to ja and e. Initially, all 
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the vertices are filled with bosons, as explained in the 
previous subsections. 

In the first order we obtain He- In the second order 
we obtain an ignorable constant, thanks to the unitarity 
of the elementary interactions. 

In the third order, as in the U(l) case, one obtains 
from the ant i- commutator contribution (see section VI A ) 
- f^^ 2 -1^ He (where J\f is the number of links, and a = 
— /i/2). From the other contribution, one obtains for each 
link 



{2\f 



{2He + {QPQ^ + QP^Q^ + Q^PQ + Q^P^Q)) 

(78) 

where the 2 comes from the bosonic creation and anni- 
hilation operators of the vertices. The MHe dependent 
terms from both contributions cancel. U sing Q PQ^ = 



e'^P.QP^Q^ = e-'^P^ {S=^, 
a renormalization to He- 



see section 



IIB2) we get 



6^'^He 



^sm 



I 



Pi + Pi 



(79) 



In the fourth order, the non-constant contributions are 
two: 

1. An "undesired" term (which is still gauge invariant, 
of course): 



H'e = i^) (0 E + Pi') (80) 



This term becomes "good" N = 2, 3, since for N = 
2, = pt2 _ I (^^Yid thus it is a constant), and 
for A/" = 3, = P"^ and P^^ = P and thus it is yet 
another renormalization of He- 

2. The anticipated plaquette terms: 

4e4- 



Hn = 



A3 



E (Qn,ig„+i,2Qi+2,i<2 + h.C.) (81) 



Thus, eventually, if we define a renormalized /i^en = 
/i ^1 — ^ sin^ (1)^ and the appropriate HE,ren = He + 

S^'^^He-! we obtain to fourth order, the desired Zat lattice 
gauge theory Hamiltonian, with some corrections: 



C. SU{N) plaquettes 

Finally, we shall describe the effective construction of 
an SU{N) gauge theory, using the appropriate elemen- 
tary interactions (for the ideal case, in which they really 
contain unitary matrices), with the loop method. Again, 
for simplicity, we describe the 2 + 1 dimensional case. We 
start with the Hamiltonian 

H = He + Hint + He (83) 

where the electric Hamiltonian of an SU{N) LGT as in 

equation ^ (again, we define M = Hint is the 

appropriate elementary interaction Q, with some fixed 
representation r of the group. We choose it to be the 
fundamental representation, and for simplicity we shall 
next drop the representation index r when referring to 
the fundamental representation. 

The spinors will be in the size of the representation - 
thus, for the fundamental representation of SU{N) we 
need N fermionic species. Besides interacting with the 
links in Hint^ they are also constrained by Hc^ which is 
similar to the previous constraining Hamiltonians. One 
such constraining Hamiltonian is 

Hc = -\ E ^t^- = ->^ E (84) 

V special v special, a 

where we define the " special" vertices to be the ones with 
both indices even: i.e., it is energetically favorable for the 
fermions to be in special vertices. As before, we intro- 
duce another set of fermions, with similar Hint, He 
constraining to other special vertices - with both the in- 
dices odd, and no interactions with the tjjs. Each plaque- 
tte contains exactly two such vertices, one of each type 
(even or odd). 

Initially, we prepare the system in the state 

Psys=p(^\{U}){{U}\ (85) 

where |{^}) ({^}| is the density matrix of some pure 
state of the bosons, and p is a mixed state, in which the 
"non-special" vertices contain no fermions at all, while 
the "special" ones are each prepared in 



^Ei«)H 



N 



(86) 



eff 



E.ren 



(82) 



(where \a) corresponds to a state of a single a fermion, tjj 
or depending on the site). 

As before, we set A to be the largest energy scale in 
the Hamiltonian, and construct an effective Hamiltonian 
for its ground sector Mq. However, as the fermions are 
here prepared in a mixed state, their "tracing out" will 
literally be tracing out, i.e. the effective Hamiltonian will 
be the result of a partial trace over the fermionic degrees 
of freedom, of the perturbative expansion, Trp {H^ffp) 



where for A' ^ oo, HE^ren He and H'^ 0, for 
N = 2 H'^ is an ignorable constant, and for A" = 3 H'^ 
is another renormalization of into He, and thus in these 
two cases the simulation is exact. 



1. Effective non-plaquette terms 

In the first order, we obtain He- Then we are left 
with an effective Hamiltonian, acting on the pure state 
of bosons. In the second order, we get once again an 
ignorable constant, due the the unitarity of the interac- 
tions. 
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The third order terms are ^e// ~ VqHiICHiICHiVo — \ {VqHiK? HiVq^V^HiVq] . Let us examine the first part. 
Consider a special vertex 1, some hnk emanating from it (in the fohowing example - negative) and the vertex on its 
other edge 2. We act with the terms constructed out of these components on a an element of a state in A^o- 

2 

^2 E (V'O (t^^)„fc(^2)j0ie2)(0ie2|i?£|0i/2) (01/21(4) f/c<i(V'i)dMi02)(rfi02| (87) 

o,o,c,a,e, J 

He does not involve any fermions, and thus (0ie2| He IO1/2) = SefHE] From the fermionic terms we get 

{^1)^ {^2), 10162) (0ie2| (^4)^ {^i)d \d1O2) {di02\ = SbeSec \a1O2) {di02\ (88) 

Considering the entire local fermionic ensemble, and tracing it out, we get -^S^c^adi and eventually the contribution 
for each negative link j^Ul^HEU^a^ and similarly, for a positive one - j^Y^UabHEUl^. 

ah 

UlHsUba = Ul {He - mE^) U^a + fiUlB^U^a = 5aa {He - mE') + fiUlE^Uba (89) 

where is the Casimir operator of the relevant link. Using the commutation relations of group elements with group 
generators (see section [ll B 3| ), we get 

E (<^- i^E - mE") + fiUlE^Uka) = MTr (A^) + NHe + {Ta)} (90) 

ab a 

where Ta is the matrix representation of Ea within the chosen representation, which for SU (N) is traceless. 

I 



For positive links, the only difference will be a change 
of sign to the last term, which is zero, so their contri- 
bution is the same. Altogether we get, neglecting the 
constant terms of Tr (A^), 4JVvj2He, where A/y is the 
number of "special" even vertices (the 4 factor is due to 
the fact, that in two spatial dimensions, there are 4 links 
neighboring each vertex). 

The other third order contribution is of the form 
—AMy^He^ and thus the entire third order contribu- 
tion (neglecting ignorable constants) is zero. This is 
all thanks to the tracelessness of the matrix representa- 
tions of SU{Nys generators. Suppose we extended our 
gauge group to U{N) - then, we would have generators 
with nonzero trace, which would lead to the Gauss's law 
terms and renormalizations, as in the U{1) case. The x 
fermions will give rise, of course, to a similar contribu- 
tion. 

The fourth order terms involving double operation of 
He, will similarly result in constants and cancelled terms, 
thanks to the tracelessness of SU{N) generators. How- 
ever, here we have one nonzero contribution, which is a 
renormalization to the electric part, 

= - E {En,k)AE^,,)^ (91) 

n,/c,o 

where for the fundamental representation of SU{N)^ 
C{N) = 1/2. 

Fourth order contributions which go back and forth 
with Hint on two links give rise to constants as well, due 
to the unitarity. 



2. The plaquette terms 

Finally, we shall construct the plaquette interactions, 
which come, as in the U{1) case, from transferring auxil- 
iary fermions around a plaquette. Here, after eliminating 
the fermionic degrees of freedom in the non-special ver- 
tices, one effectively gets, for example, terms like 

- ^ E W)a (f^l)a6 (f^2),e {U^) (f/])^^ (^„), |e) (e| 

abcde 

(92) 

(with the conventions of figure [sj with v the vertex la- 
beled there by n). By tracing out the fermions, we get 
Sae, which "closes" the plaquette and introduces the de- 
sired group-trace. This yields the plaquette Hamiltonian 

9 4 

^s = -^ E (Tr(c/iC/2C/]C/i) +/I.C.) (93) 

plaquettes 

(the factor 2 is due to the two types of auxiliary fermions) 
and eventually we get effectively, up to fourth order, the 
SU (N) pure-gauge Kogut-Susskind Hamiltonian, 

Heff=HE^SHE^HB (94) 

without any corrections {SHe is just a renormalization)! 
Moreover, the leading order of e is the fourth one, and 
thus here it is even sufficient to demand <C A^. 
One should also note, that by introducing other set of 
fermions, interacting with the bosons with a similar Hint 
but non-constrained, one can introduce dynamic fermions 
to the system. 
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VII. QUANTUM SIMULATION OF 2 + 1-D 
GAUGE THEORIES 



Putting together the elementary link interactions of 
section |IV[ and the loop method for the plaquette in- 
teractions of section |VI[ one could construct quantum 
simulations of lattice gauge theories in 2 + 1 dimensions. 
However, in several cases one would have to face "real" 
conditions, instead of the "ideal" conditions of the previ- 
ous section, such of the use of finite Hilbert spaces instead 
of infinite ones. 

In that sense, the quantum simulation of Zat is the 
most accurate: in this theory, the local hilbert spaces are 
already of finite dimension, and thus nothing should be 
added to the discussion of these models in the previous 
section. 

On the other hand, quantum simulation of the contin- 
uous gauge theories - U{1) and SU{N) - require some 
truncation of the Hilbert space (This was already men- 
tioned in the case of the Schwinger model, in section 
VA). In the case of cQED, the simulation requires the 
use of BECs. In this section we show that it should still 
work with a finite number of bosons per link, derive the 
conditions for that, and give numerical evidence. Finally 
we shall comment on the SU (2) 2 -f 1-d simulation. 



2 + 1-d simulation of cQED using a finite 
number of bosons 



We have shown in section VIA how to obtain, ef- 



fectively, the plaquette interactions using the already 
gauge invariant elementary interactions. However, in the 
derivation we have used unitary matrices in the elemen- 
tary interactions {Nq oo), whereas the in the "real" 
scenario Nq = £ is finite, even if it's large, and thus the 
interactions contain angular momentum ladder operators 
( [38| ) which are nonunitary. Here we shall describe how 
to handle these effects in order to achieve an accurate 
simulation despite the non-unitarity. 

We still include the x fermions, with the same con- 
straining part fnl, but now Hint takes the form 



Hint — 



A £(£ + 1) 



T 



± 



£{£^1) £{£^1) 



(where Lj^L_ is for positive links, and vice-versa for the 
negative ones). This contribution becomes constant as 
^ ^ oo, because then, any eigenvalue of Lz satisfies m <C 
£. In the case of a finite ^, one can neglect the first 
constant, and after summing on all the links, obtain a 
renormalization factor for the electric Hamiltonian from 
the term, and a Gauss's law (which is an ignorable 
constant) from the linear Lz part which has the correct 
signs. Thus, second order leads to a renormalization of 
He (and this is why the xs are important). 



2. Third and fourth order contributions 

Here, after cancelling the equal part in the two possible 
contributions (see the previous section), we are left with 



(l^2L^ ^^^^ 



(97) 



where the choices of signs for positive/negative links are 
as in the second order. 

Again, in the ideal limit, j^xr\ — ^ 1? ihen we 

get a constant + Gauss's laws (which are constants too). 
However, for a finite ^, the contributions are again non 
constant: besides the constants, we get linear terms in 
Lz, which correspond to Gauss's law, and a third order 
renormalization of He, but also L^ terms. If we work 
with £ = No = 1^ L^ = Lz and thanks to the correct 
signs get Gauss's law again. However, this is not a very 
interesting and obviously not the general case, and thus 
we wish to find some way to deal with these extra terms. 
We shall first focus on the fourth order contributions, 
and then conclude what to do about the "problematic" 
terms from both orders. 

In the fourth order there are three types of contribu- 
tions. First, the plaquette terms, which now take the 
form 



Hn = 



2e^ 



(^+,n,li^+,n+i,2^-,n+2,l^-n,2 + h.c}j (98) 



h.C.) 



(95) 



L First and second order contributions 

In the first order, we get He as before. In the second 
order, we get the same type of contributions, but now 
they will not be constant anymore. What we have now 
are contributions of the form 



(96) 



The second type of terms involve operations with Hint 
only, and they include, after the reduction of constants, 
products of Lz and L^ of neighboring links intersecting 
in odd vertices. This term is in the same order of the 
plaquette terms - O (^js^. These terms, which we call 

i^^, are, indeed, unwanted terms, but can still be toler- 
ated, as they are gauge invariant and not stronger than 
the plaquette interactions. 

The ^ third terms take the form 
~ xHii+i) {^z^L±}{Lz,L^} on positive/negative 
links. These terms will renormalize He in the ideal 
limit, but for a finite £ will include, besides the constants, 
Lz and L^ terms, which are "treatable", also Ll and 
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FIG. 14. A plot of d (equation ( |1Q4| )) for several values of One can see that d grows as decrease. However, in this very 
same region of parameters, the credibility of it is low, as explained in the text. 



terms (again, with the correct signs to contribute to 
Gauss's law and renormalize He for 1 = 1). In order 
to eliminate the effect of these terms, along with the 

2 

undesired third order terms, we define /i = where 
/3 < 1 is a dimensionless parameter. Then, the undesired 

third order terms are O ^-^-^ - like the plaquettes, 

and thus can be tolerated, at least for states which 
are superpositions of m <C ^ mostly. The fourth order 
undesired terms become effectively sixth order terms - 

O - and thus can be safely neglected. 

Thus, for a finite A^o = if the parameters are tuned 
correctly, one gets effectively, up to constant, the Hamil- 
tonian 



H^ = He^Hb^O 



A3 



with 



He 



(99) 



(100) 



n,/c 



which is expected to give rise to the same dynamics as 
the Kogut-Susskind Hamiltonian H = He-\-Hb, at least 

4 

for a regime in which is small enough: i.e., not in the 
extreme weak limit, but apparently not only in the strong 



limit, but also in a regime where He^ Hb are of the same 
order of magnitude. The two Hamiltonians have to differ 
by a constant, at least for states which are superpositions 
mostly of m <C ^. 

In order to see that, we shall consider some numerical 
results. But before that, let us point out two important 
issues about the Hamiltonian parameters. Although the 

2 

relation /i = /3y might seem to imply that He will al- 
ways be stronger than Hb^ because of the powers of e and 
A, one should note that thanks to /3, which can be chosen 
small, one can still go to weaker regimes. Moreover, as 
in the ideal case. Hint only contributes to the effective 
Hamiltonian in even repetitions, and thus, again, it is 
sufficient to demand <C A^ rather than e <C A. 



B. 2 + 1-d cQED simulation - A Numerical proof of 
principle 

In order to see whether the "real" case simulation still 
yields valuable results, one has to check the effects of the 

O (^fs^ in (99). This can be done if one compares the 

spectrum of the'^ desired" Hamiltonian, 



Hi = He + Hb 



(101) 
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FIG. 15. Plots of the spectra of H (upper panel) and H (lower panel) for small values of I and several values of x, compared 
to the fairly reasonable approximation of the non-truncated Kogut-Susskind model using I = 100. Note that the qualitative 
correspondence of the two spectra arises as x,^ increase. 



with the spectrum of the lower energy sector (i.e., the 
states fulfilling the constraint) of the fundamental Hamil- 
tonian, 



H = Hf 



He 



(102) 



The accuracy of the simulation will be deduced from 
the observation of a constant energy shift between the 
spectra of the two Hamiltonians. This will mean that 
both of them give rise to the same dynamics. Such results 
will trivially take place in the strong coupling limit, where 
e ^ 0. However we wish to check what happens in other 
coupling regimes. 

We have run numerical simulations of a single plaque- 
tte, using the parameters A = 10, e = 0.1. We changed 
the values of ^, /3, using the convenient dimensionless ra- 
tio between the scales of He and Hb^ 



X — 



13 



1 



2 VA £{£ + !) 



A2 



(103) 



parametrizing the coupling strength. We have chosen the 
subspace with no (real) static charges. 



We have calculated the spectra {E} of i^'s ground sec- 
tor, as well as |^| of ^ for several values of £ and x and 
calculated the value of 



std 



mean [E — E 



(104) 



as one can see in figure [T4j this value is fairly small, get- 
ting larger, as anticipated, for smaller values of the sim- 
ulation parameters. However, on the other hand, if one 
wishes to evaluate the quality of simulations, it should be 
noted that in order to get an accurate simulation (com- 
pared to the original Kogut-Susskind model) in the weak 
coupling limit (small x) one has to pick a sufficiently 
large value of £ [13]. If only the ground state is of in- 
terest, lower -^s are possible. In order to achieve a better 
accuracy for excited states, one must increase £. This can 



be seen in figure 15 



Furthermore, as x gets smaller, the energy scales of 
the two Hamiltonians H^ H separate (as can be seen in 
figure 15). Thus, the difference between them can still 



be treated approximately as a constant, due to the dif- 
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FIG. 16. Plots of the spectra of H,H, enabling a qualitative impression on the constant difference between them, for several 
valus of the simulation parameters x/. 



ference of scales, but yet one should note this qualitative 
subtlety as well, treating the left region of figure [M] more 
carefully. On the other hand, for x > 1, using sufficiently 
large ^s, which tend to get smaller as x is increased (see 
figure [l6|, the scales of two Hamiltonians correspond and 
the d value is credible. As one can conclude from that. 



a simulation which is accurate both comparing to the 
Kogut-Susskind model and with fairly constant effective 
Hamiltonian contributions is possible for x > 1. 
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C. 2 + 1-d SU{N) simulation 
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In the case of non-abelian theories, such as SU(2)^ re- 
lying on [17 for the elementary interactions (as described 
one could simulate a version of SU{2) 



IVC) 



m section 

gauge theory in which the bosonic operators in the ele- 
mentary interactions are not unitary (63). 

In the 1 + 1-d case we could establish a good approxi- 
mation for the unitary interactions, since the elementary 
interactions only change by ±1, or j by ±1/2 on single 
links. However, here, since there are plaquette interac- 
tions, we have such processes on four links at once. One 
can also introduce dynamic fermions, and altogether we 
lose the fifth-order accuracy we had in the 1 + 1-d case. 
Moreover, the effective Hamiltonian series may include, 
in this case, more terms (as we had in the "real" cQED 
case comparing to the ideal unitary case). This will de- 
pend on these TVs as well. On the other hand, note that 
all these "corrections" are nevertheless gauge invariant - 
gauge invar iance is not ruined! 

Therefore we conclude that due to the \/N + 1 opera- 
tors, which would lead, already in first order in (p3|, 
to terms with different amplitudes than in regular SU(2) 
theory, the [17 realization of SU(2) elementary interac- 
tions allows only for simulation of SU{2) gauge theories 
near the strong coupling limit. However, in this regime, 
it should reproduce the same physics, qualitatively. 



VIII. SUMMARY AND CONCLUSIONS 

In this paper we have presented a new method for 
the simulation of lattice gauge theories in high energy 
physics, using ultracold atoms in optical lattices. Unlike 
previous proposals, where the gauge invar iance was ob- 
tained effectively, here the gauge invariance is fundamen- 
tally in the atomic Hamiltonian, without using perturba- 
tion theory. This is done by utlizing the conservation of 
hyperfine angular momentum, introducing the gauge in- 
variant elementary matter-gauge field interactions along 
links. 

At a second stage, the nonlinear self-interactions of 
the gauge fields along plauqettes are constructed effec- 
tively, using the introduced loop method, from the al- 
ready gauge invariant Hamiltonian with the elementary 
interactions. Although the derivation uses fourth order 
perturbation theory, it is effectively to second order, due 
to the vanishing of the relevant odd orders. 

Along the simulation proposals in these paper, was the 
suggestion to simulate compact QED in 1 + 1 dimensions 
- the Schwinger model. This simulation scheme utilizes in 
this case a small number of atomic species and does not 
involve perturbative and effective methods at all. Thus 
it could serve as the first step towards the realization of 
quantum simulations of gauge theories and high energy 
physics models. 
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APPENDIX: FURTHER DETAILS ON THE 
LATTICE GAUGE THEORIES DISCUSSED IN 
THE PAPER 

1. Confinenemt in cQED 

In this appendix we briefly review confinement in 
cQED^ whose Hamiltonian was discussed in section 

ittbti 

Let us consider first the states of the pure gauge the- 
ory, i.e. with no dynamical matter. In this case only 
static charges are possible, and the Hamiltonian is the 
abelian Kogut-Susskind Hamiltonian (14). Suppose we 



use the local flux basis - i.e., the basis consisting of prod- 
ucts of flux states over the links of the entire lattice. In 
the strong coupling limits where ^ oo, this basis is a 
good choice, as its elements are eigenstates of He- Note, 
however, that these eigenstates are divided into several 
sectors, depending on the static charges - eigenvalues of 
Gn, as these are constants of motion, setting the Gauss's 
law (13). Let us discuss this limit, taking Hb as a per- 



turbation to He- 

Consider, in the strong limit, the case of a single charge 
Qn = 1 in some vertex n. In order to respect Gauss's 
law, an infinite string of flux must be introduced to the 
system. In fact, the (degenerate) zeroth order ground 
state of this charge sector is already of infinite energy: 
an infinite string of flux ±1 must emanate from n. This 
is an evidence of charge confinement in the strong limit: 
a single charge "costs" infinite energy. Thus we must 
consider bound states of at least two charges. Let us 
have a look on the state of two such charges, say = 1 
and Qn+i?k ~ There, the zeroth order eigenstate 

consists of a string of flux 1, connecting the charges, i.e. 
the state 



R-l 



(105) 



where V^j^, i/j^^ create positive/negative static charges in 
the appropriate vertices and \vac) is the zeroth order 
chargeless ground state, with flux zero everywhere. Note 
that this state is gauge invariant (just check t he a belian 
version of the local gauge transformation laws (SjS ). This 
state describes an electric flux tube, as expected within 
a confining phase. The energy of such a state, to zeroth 
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order, is just the eigenvalue of He'- 



Hp 



2 



(106) 



- the static energy of this particle-antiparticle bound pair 
(a "meson") is proportional to R (up to higher order 
corrections, O • This manifests the hnear potential 
law, expected in a confining phase: V {R) oc R^ unlike 
the well-known Coulomb phase behavior, V (R) (x R~^ . 

The weak limit g'^ 0, however, is different. There, 
Hb is stronger, and these fiux states are no longer eigen- 
states of the system, even not approximate ones. In 3 + 1 
dimensions, there is a phase transition to a Coulomb 
phase [26 . In 2 + 1 dimensions, the theory confines 
for all values of the coupling constant, and it is a non- 
perturbative effect [24l I33H35] . This is also the case for 
1 + 1 dimensions, where the theory is exactly solvable. 
When temperature enters the game as well, there is a 
phase transition to a Coulomb phase for small ^s also in 
the case of 2 + 1 dimensions [69l [70] . 

One may wish, as we do here for simulation purposes, 
to consider the case of a finite, truncated, local Hilbert 
space. In this case, the strong limit behavior is unaf- 
fected. For the weak limit, only when it is confining one 
can approximate the behavior of the real theory - as in 
the confining phase only a few number of low fiux states is 
needed. The approximation becomes better as the trun- 
cation value ^ is stronger, but also for small values of ^ a 
reasonable approximation holds [13]. This is true when 
one considers the lattice theory; In terms of the contin- 
uum limit, if the truncation is done properly one can 
regain the right continuum limit - see, for example, the 
link models 



2. Dynamic fermions in cQED 

Including dynamic fermions in a lattice gauge theory 
is a bit problematic, due to the doubling of fermions in 
the continuum limit. We focus, in this paper, and thus 
here, in the method of staggered fermions [36j. There, 
the continuum spinors decompose into two neighboring 
vertices of the lattice, one containing particles and the 
other - antiparticles. 2-component continuum spinors 
(which are sufficient for 1 + 1 dimensional theories, and 
also for 2 + 1 dimensional ones, if one is not intersted in 
the chiral anomaly) are thus formed of two neighboring 
vertices, containing up to a single fermion each. This is 
why the Hm part of the Hamiltonian (35) contains the 



masses with alternating signs: the odd vertices represent 
the antiparticles, and they are filled in the ground state 
- forming the Dirac sea. Moreover, the dynamic charge 
Qn, present in Gauss's law ^ is defined as 

gn=ViV'n-^(l-(-ir+"^) (107) 

(for simplicity this is a 2 -hi dimensional definition, which 
can be generalized to other dimensions, keeping the ver- 



tices' parity properties). This is in with the particle/anti- 
particle picture: even vertices, representing particles, can 
be occupied by no fermions - corresponding to no mass 
and no charge - or a single fermion, corresponding to a 
particle with mass M and charge +1; Occupied odd ver- 
tices have no mass (relative to —M) and no charge, and 
correspond to vacant vertices in the HEP picture, while 
empty vertices have mass M (relative to —M) and charge 
— 1, representing an anti-particle. 

The fermionic interactions are obtained with the gauge 
invariant i^int, whose abelian form is (15). However, one 



should note that in order to get the Dirac equation in 
the continuum limit, some phases should introduced to 
this Hamiltonian, and they can be achieved using canon- 
ical transformations (which are not required to be imple- 
mented experimentally) on the fermions. For example, in 
a 1 + 1 dimensional system, one can use the transforma- 
tion in [64 . The phase prescription for 3 + 1 dimensions, 
yielding four- component spinors, is given in [36 . 



3. Confinement in Zat gauge theories 

Confinement for these theories is especially interesting, 
as it was argued [71 that confinement of quarks in QCD 
{SU{3)) has to do with the group's center, which is the 
group Z3. Thus the phase structure of Z^v is interesting 
j37l \72\ [73] . The theory confines in the strong coupling 
limit {jii 00). Both in 2 + 1 and in 3 + 1 dimensions 
it in not the only phase. In 3 + 1 dimensions, due to the 
self-dual nature of the theory, there is a phase transition 
from electric confinement for large coupling to a magnetic 
confinement in the small coupling regime, for N < Nc 
{Nc ~ 6). For N > Nc there is a third phase, with no 
confinement at all. In 2 + 1 dimensions there is a phase 
transition and the theory does not confine in the weak 
limit [37] • However, for A/" ^ 00 the theory shows the 
phase transition at ^ = [73 , in accordance with the 
single confining phase structure of U{1) (the ^ 00 
limit of Zat). 



4. SU{N) gauge theories - dynamic fermions and 
confinement 



In order to introduce dynamical fermions, we choose 
to utilize the method of staggered fermions again. We 
choose, as usual, to suppress fiavor and spin indices, and 
identify the spinor components only by group indices, 
as in equation ([T]). With staggered fermions the masses 
are positive/negative alternately, similarly to (35). This 
time, the non-abelian charge is defined as 



{Qn)a = E (Ta),, (^„)e = (^i)^ K)6c (V-n), 



be 



(108) 

(where the second equality holds for the fundamental rep- 
resentation). This is a fermionic Schwinger representa- 
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tion, and thus a full vertex has a zero charge. Thus one 
can establish a Dirac sea, similarly the abelian one (ap- 
pendix 1), where vacancies in odd vertices correspond 
to anti-particles, measuring the masses there relative to 
— M per fermion. 

One can also consider static charges (i.e. with infi- 
nite mass). In particular, one would like to examine 
the confinement of static fundamental charges. This was 
initially done by Wilson [23] and later other authors in 
the Euclidean approach [24l [33], but we shall concen- 
trate on the Hamiltonian picture. Consider the strong 
limit of the theory, ^ 1, in which the electric part 
is much stronger than the magnetic part, which may 
be thus considered as a perturbation (as in the abelian 
case). In this context, the ground state of two charges in 



n, n + i?x, is the non-abelian analogue to (105): Thus, 



the links along the straight line connecting the charges 
are excited to j = ^, contributing zeroth order energy 

\j (j + 1) = \\i while the other links are in the sin- 
glet j = 0, contributing no energy at all in zeroth order. 
Therefore the energy of this state, to zeroth order, is 
E {(/)) = \ \R (as there are R excited links), and we get 
confining behaviour, V (R) oc R again: the static energy 
(potential) between two static "quarks" is proportional 

to their distance. The string tension, E/R, is |^ in 
this case [27 . Other coupling regimes will not be treated 
here. We shall only comment that as the coupling de- 
creases within a confining phase, the fiux tube broadens 
and 1^) is no longer an eigenstate of the Hamiltonian, in 
accordance with the growth of perturbative corrections, 
until the breakdown of perturbation theory. 
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